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Abstract. — We consider the following coupled cubic Schrodinger equations 
00 ( id t u + d 2 u = e 2 \v\ 2 u, (i^jslxS 1 , 

^ idtv + d 2 v — e 2 \u\ 2 v . 

We prove that there exists a beating effect, i.e. an energy exchange between different modes. This 
Oh construction may be transported to the linear time-dependent Schrodinger equation: we build solutions 

such that their Sobolev norms grow logarithmically. All these results are stated for large but finite times. 

Resume. — Nous considerons le systeme d'equations de Schrodinger couplees 

f id t u + d 2 u = e 2 \v\ 2 u, (t,x)€RxS 1 , 

> 

o 

oo 
\o 
m 

00 

o 

(N 



id t v + d 2 v = e 2 \u\ 2 v. 

Nous montrons l'existence d'un effet de battement, c'est-a-dire un echange d'energie entre des modes 
differents. Cette construction peut etre transposee pour l'equation de Schrodinger lineaire non autonome, 
ce qui permet de construire des solutions dont les normes de Sobolev croissent logarithmiquement (inflation 
de normes). Tous ces resultats sont etablis pour des temps grands mais finis. 



1. Introduction 



1.1. General introduction. — Denote by S 1 = IR/27rZ the circle, and let e > be a small 
parameter. In this paper we are concerned with the following cubic coupled non linear Schrodinger 
equations 

' idtu + dlu = e 2 \v\ 2 u, (^)elxS 1 , 
(1.1) < id t v + d 2 x v = e 2 \u\ 2 v, 

k u(0,x) = u Q (x), v(0,x) = v (x) . 
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We exhibit some solutions of this system which stay close to solutions of a finite dimensional 
nonlinear system for long times. We stress out that these solutions are not obtained by perturbations 
of the associated linear system. Thanks to the nonlinearity, we may produce a beating effect, i.e. a 
transfer of energy between two different modes, something which is not possible in the linear case. 
The solutions of the initial system are then found thanks to a resonant Birkhoff normal form and 
approximation arguments, and they enjoy the same beating properties as those of the reduced system. 
This phenomenon heavily relies on the presence of resonances. Actually, Bambusi and Grebert [2\ 



showed that, in the non-resonant setting (e.g. adding a typical potential in each equation of (1.1)), 
the dynamics stays close to linear for long times (see the introduction of [5]). 

This new example leans on a principle that was already used in [5] and [5]: we make it explicit in 
Section 2, where the conditions for applying our method on different resonant Hamiltonian PDEs are 
enumerated. 

The control of Sobolev norms in Hamiltonian PDEs has a long story, both in the nonlinear and the 
linear time-dependant setting. Concerning nonlinear equations, one of the most outstanding results is 
due to [5 J for the cubic 2-dimensional NLS, recently completed by [10] . where there is a construction 
of specific solutions which exhibit a polynomial growth of Sobolev norms for large finite times. 

In the linear setting, Bourgain |3] proves a polynomial bound of the Sobolev norm of the solution 
u of 

id t u + d%u + V(t,x)u = 0, (i.ijElxS 1 , 

where V(t, x) is a bounded (real analytic) potential. Moreover, when the potential is quasi-periodic 
in time he obtains in [4j a logarithmic bound. This last result has been enhanced by Delort [6] and 
Wang |12j . who gets a logarithmic bound for bounded potentials. As a by-product of our work, we 
recover a result of Bourgain [4], who showed that these logarithmic bounds are optimal in the case 
of analytic potentials (see Section [4]) . Note that it is possible to obtain a growth of higher order 
(but still logarithmic) when considering potentials in Gevrey classes (as in Fang-Zhang |7J), or even a 
sub-polynomial growth in the case of C°° potentials. 



1.2. Beating effect in the system (1.1). — Our first result concerns the dynamics of (1.1). 



Theorem 1.1. — For all < 7 < 1/2, there exist < T 7 < C|ln7|, a 2T 7 — periodic function 
Kj : M 1 — >]0, 1[ which satisfies K~(Q) = 7 and Ky(Ty) = 1 — 7, and there exists Eq > so that if 



p,q 6Z and if < e < min(eo,7 )> there exists a solution to (1.1) satisfying for all \t\ < e 



--3 



u(t,x) = Upi^jP* + u q (t)e i i x + e 1 / 2 r u (t,x), 
v(t,x) = v p (t)e ipx + v q {t)e iqx + e l / 2 r v {t,x) , 



with K(*)l = \ v p{t)\ = K 7 (e 2 t 

\u p (t)\ 2 = K(*)l 2 = l-K^e- 



and where r u and r v are: 

- smooth in time and analytic in space on [— e -3 ,e~ 3 ] x S . 

- for r = r u , r v the Fourier coefficients rj(t) of r(t) satisfy for some p > 



sup |rj(i)| < Ce 

|t|<e-» 



-P\j\ 
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uniformly in e > and p,q G Z. 

This statement shows an exchange of energy between the modes p and q: the mode u q and v p grow 
from 7 to 1 — 7 in time t = e~ 2 T^. Considering the larger time scales e~ 3 , we obtain a periodic 
phenomenon which we will call beating effect. 

Of course the solutions satisfy the three conservation laws: the mass, the momentum and the energy 
are constant quantities. 

• Conservation of the mass: J \u\ 2 and J \v\ 2 

(1.3) |ug| 2 + \u p \ 2 = est, \v q \ 2 + \v p \ 2 = est. 

• Conservation of the momentum: Im J ud x u + Im J vd x v 

(1.4) q\Uq\ 2 +P\Up\ 2 + q\v q \ 2 +p\v p \ 2 = CSt. 

• Conservation of the energy: J \d x u\ 2 + |<9 x -u| 2 + e 2 J 

fi c\ 2i |2 I 2 1 1 2 i 2i |2 i 2 1 1 2 j. 

(1.5) q \u q \ +p \Up\ +q \v q \ +p \v p \ = est. 

On the other hand, the solutions given by Theorem 1 1 . 1 1 satisfy for < t < and s > 

(1.6) \\u(t,-)f H s = (q 2s - P 2s )K 7 (eH)+ P 2s + 0(e). 

In particular, this norm does not remain constant in time, which is a true nonlinear effect. However, 
the sum \\u(t, -)\\% s + \\v{t,')\\%a remains almost constant and thus ( |1.6[ ) cannot be interpreted as a 
norm inflation. Nevertheless this effect will be used in the linear case (cf. Theorem 1.3). 

Remark 1.2. — For the defocusing- focusing system 

' idtu + dlu = e 2 \v\ 2 u, (^)elxS 1 , 
< idtv + d 2 v = —e 2 \u\ 2 v, 
k u(0, x) = u (x), v(0, x) = v (x), 



, 1 2 1 \2 
\U\ \v\ 



(1.7) 



one can show a beating effect only for the case q = —p (see also Remark 3.3). 

1.3. Growth of Sobolev norms in linear Schrodinger equations. — Theorem |1.1| allows us 
to build real time-dependent potentials V(t, x) for the following linear Schrodinger equation 

(1.8) idtu + d 2 x u + V(t, x)u = 0, (f.ijelxS 1 . 



Namely, in (1.1) we consider the solution v as a given function and we set V(t, x) = — e 2 \v (t, x)\ 2 . 
For a > 1 we define the Gevrey class ^(S 1 ) as the set of functions / G C°°(S 1 ) satisfying, for some 
A > and C > 0: 

sup \f {n \x)\ < CA n (n\) a , VraeN. 

xes 1 

In the periodic setting, an equivalent formulation is available (see |llj ): a function / G (^^(S 1 ) is in 
^(S 1 ) if, for some K > and B > 0, we have for any j G Z, 

(1.9) \f 3 \<Ke- B \^\ 
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where (fj)jez denote the Fourier coefficients of /. We then define a semi-norm as the best 

constant K in ( |1.9| ) (see for more details). 

The beating phenomenon then leads to the growth of Sobolev norms (for finite but arbitrary large 
times) for some solutions of this equation. Obviously, since V is a real potential, the L 2 norm of any 



solution of (1.8) is constant. However, we are able to prove 



Theorem 1.3. — Fix s > and a > 1. There exist a sequence of real potentials V q (t,x), a sequence 
of initial conditions (uq) and a sequence of times T q — > +oo as q — > +oo such that 

- The potentials are smooth in time, real analytic in space and uniformly bounded in Gevrey classes: 

V?a, ViG [0,T q ], \\V q (t, .)\\g a < C a , 

- \\ u o\\h° = 1, 

- The corresponding solutions to the Cauchy problem u g (t, .) are real analytic in space fort G [0, T q ], 

- There exists a constant C a ,s depending only on a and s such that 

\\u q (T q )\\ H s > C a , s (lnT q ) sa . 

This can be compared to the result obtained in the analytic case (both in time and space variables) 
by Wang |12j . who proves that if the potential V is real analytic in (t,x) in a band D = (M + ip) 2 , 
real and bounded in M 2 , then, given any s > 0, there exists C s such that 

\\u{t)\\ H s <C s {\og{\t\+2)Y s \\u \\ H s, 

where k is a constant independent of s,uo. 

Bourgain also showed in [4] that < C(lnt) c ' s ||«o||j^s when V is analytic and quasi-periodic 

in time. This has been extended by D. Fang- Q. Zhang [7]. If V is C°° but not supposed to be 
quasi-periodic in time, Bourgain [3] proves the bound ||u||_f/ s < Crp \ \ "Uoll-ff 3 ; f° r au V > 0- See also 
the nice generalisation by J.-M. Delort [6] . 

1.4. Plan of the paper. — We describe in Section [2] the normal form method used to extract 
from the infinite dimensional Hamiltonian system a nonlinear finite dimensional and integrable system 
which will drive our solutions. The solutions of this small system are then studied in Section |3| The 



proof of Theorem 1.1 then relies on the control of the other terms in the initial Hamiltonian system, 



that is the topic of Section HI Finally in Section [5] we prove Theorem 1.3 



2. The normal form 

2.1. Principle of the result. — In this section, we formalise the principle already used in [9] and 
[8], in order to follow for arbitrary long times solutions of an integrable model equation. The system 
has to be Hamiltonian: let H be the Hamilton function describing its dynamics on some Hilbert phase 
space. We assume that H is smooth in a neighbourhood of the origin, and that its Taylor expansion 
is given by 

h = n + z; es + z; r + r p+1 , 

where 
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N is a homogeneous polynomial of order 2, coming from the linear part of the system. It usually 
gathers the linear actions, i.e. the first integrals of the linearized system at the origin, which may 
be easily written in action-angle coordinates thanks to a Fourier transform for instance. We take 
the following form for clarity: 

where the \j are the eigenvalues of the linearization at of the system. We suppose that A-,- grows 
polynomially with j: Xj ~ \j\ r , with r > 1. 

For p an even integer, Z p es + Z p r is the next nonzero term in the Taylor expansion of H. It is a 
homogeneous polynomial of degree p. We distinguish between resonant and nonresonant terms in 
the sense of Birkhoff normal forms: a monomial M of degree p is called resonant if it commutes 
with N, i.e. { M, iV } = 0. On the contrary case, a nonresonant monomial may be removed by 



one step of Birkhoff normal form (see Proposition 2.1): we suppose that a Birkhoff normal form 
is available for the system in a ball B centred at the origin. 

- Rp+i is an analytic Hamiltonian which vanishes at the origin up to order p + 1. 

In order to observe some beating effect, we have to focus on the resonant part. Suppose that we may 
decompose 

N + z; es = H n + N ext + Z p ,i + z pa + z p>2 , 

where 

- H u , defining the reduced Hamiltonian system, depends only on finitely many variables (indexed 
by j £ A), called the internal modes, 

- N ext contains all the monomials of N depending on the external modes, i.e. the variables indexed 
by j A, 

- Z p i gathers monomials involving exactly one external mode, 

- Z P) 2 gathers monomials involving exactly two external modes, 

- Z p > 2 gathers monomials involving at least three external modes. 

We may now write the principle already used in [9j and [5], put in light again in this paper. This 
brings together the assumptions needed to exhibit beating phenomena for Hamiltonian PDEs using 
our method. 

Principle. — If the following assumptions are fulfilled: 

- H u defines a completely integrable Hamiltonian system, 

~ 1 1 — y (q(t),p(t)) is a solution of the reduced system satisfying that for every t, (q(t),p(t)) stays in 
the ball B, 

- Z Pj i = 0, i.e. resonances cannot light on one single outer mode, 

- { N ext ,Z Pi 2 } = ; i.e. Z Pj 2 does not affect the external modes. 

- There exists a strictly convex combination of the (Ij)je.A denoted by I such that {l,H n } = 0. 
Then there exists solutions of the system governed by H which follow (q(t),p(t)) for long times, i.e. 
their projection on the reduced phase space stay close to (q(t),p(t)) and the difference between the 
solution and its projection stays small, for long times. 

The beating effect is then obtained when we are able to construct a periodic solution t h-> (q(t),p(t)) 
of the reduced system. Note that the 2D cubic NLS equation enters in this setting when considering 
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"small squares" of indices, e.g. I = {(0,0), (1,0), (0,1), (1,1)} in the Fourier modes decomposition. We 
do not write the details. 



2.2. Hamiltonian formulation and Birkhoff normal form. — To apply a normal form proce- 
dure it is convenient to transform the original system where the nonlinear term is small into a system 
where the solutions are small. Namely by an obvious change of variable, (1.1) is equivalent to the 
system 

' id t u + d 2 x u = \v\ 2 u 1 (t.ijeKxS 1 , 
< idfV + d 2 v = \u\ 2 v, 

k u(0, x) = £Uq(x), v(0, x) = £Vq(x) . 



(2.1) 

Denote by 



H 



I 



d T u\ 2 



\d x v\ z + I \u\ 2 \v\ 2 , 



the Hamiltonian of (2.1) with the symplectic structure du A du + dv A dv . In other words, (2.1) is 
equivalent to 



(2.2) 



u 



,6H_ 
5u ' 
.5H_ 

8v ' 



u = i 



V = I 



Let us expand u, u, v and v in Fourier modes: 



u 



,6H_ 

5u ' 
5H_ 

5v 



a 9 e 



W = E 

(x) = J2^ ijx , v(x) = J2Pje 



-IJX 



-IJX 



We define 



P(a,(3) 



/ \u(x)\ 2 \v(x)\ 2 dx = V a jx <Xj 2 Pt. x P t , 

.1Z,0. 



M(j,£)=o 



where Ai(j,£) = j\ — j2 + i\ — (-2 denotes the momentum of the multi-index (a, (3) 6 Z 4 . 
In this Fourier setting the equation (2.2) reads as an infinite Hamiltonian system 



(2.3) 



iotj = j ctj + 



dP 


dH 


daj 




dP 


dH 


w 3 


~~w 3 



dP dH 



-la 



dctj dctj ' 



, 2 - + dP = dH 

df3 3 5/3/ 



For p > 0, we consider the following phase space 



T p = { (a,/3) G (f(Z))\s.t. \\(a,P)\\ P ■= E^N + l&D < 00 
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(2.2) is equivalent to 



which we endow with the canonical symplectic structure —i ^^(day A dotj + d(3j A df3j). According to 

j 

this structure, the Poisson bracket between two functions / and g of (a, a, (3, 0) is defined by 
If \ — ■ V r d f dg df dg , df dg df dg > 

It is convenient to work in the symplectic polar coordinates (ay = y^e 4 ^ , ay = \/Tje~ %e \ j3j = 
/ ~Jje lipj , f3j = -\J~Jje~ lVj ) - eZ - Since we have da A da = idO A dl and d/3 A d/3 = id(/2 A dJ the system 

• dH ■ dH 

dH ■ dH 

We denote by B p (r) the ball of radius r centred at the origin in T p , and introduce the resonant set 
n = { (jl,h,h,h) G Z 4 s.t. | ji -j 2 + e 1 -£ 2 = and - j 2 2 + i\ - l\ = } . 

Proposition 2.1. — There exists a canonical change of variable r from B p {e) into B p {2e) with e 
small enough such that 

(2.4) H:=Ho T = N + Zi + Rq, 

where 

(i) N is the term N{I) = /^j 2 (Ij +Jj)- 
(ii) Z± is the homogeneous polynomial of degree 4 



Z 4 = ^a h a j2 l3 h P 



ti- 



ll 



In particular, Z4 is made of resonant monomials: it satisfies { Z4, N } =0. 

(Hi) i?6 is the remainder of order 6, i.e. a Hamiltonian satisfying 
\\X Re (z)\\ p <C\\z\\ 5 p forz = (a,a,(3j)GB p (e). 

(iv) t is close to the identity: there exist a constant C p such that \\t(z) — z\\ p < Cy,||;z[|? for all 
z € B p (e). 

The proof is similar to the proof of Proposition 2.1 in [15]: we essentially use that if A4(j,£) = 
and (j, I) U then \ j\ - j\ + 1\ - l\\ > 1, i.e. there is no small divisors involved. 
For the construction of a more general Birkhoff Normal Form see [2]. By abuse of notation, in the 
proposition and in the sequel, the new variables (a',/3') = r _1 (a,/3) are still denoted by (a,/3). 
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2.3. Description of the resonant normal form. — In this subsection we study the resonant 
part of the normal form given by Proposition |2.1[ Denote by 



I= J2\ a n\ 2 , J = ^2\Pn\ 2 , S=^2a n /3 n . 

n&L n£Z neZ 

Proposition 2.2. — The polynomial Z4 reads: 

Z 4 = IJ+\S\ 2 -Y,\<*n\ 2 \f3n\ 2 . 

Proof. — By an elementary computation, we know that (ji, j'2,^1,^2) £ 1Z iff {j'1,^1} = {.72,^2} and 
the result follows. □ 



3. The reduced model 



We want to describe the dynamics of the Hamiltonian system obtained by reducing (2.3) to the 
space 

J{P, q) ■= { (a, 0) G ?p I ay = ay = =p j = when j ^p,q}, 
and we denote by H the reduced Hamiltonian, i.e. 

H = H\ 



\J(p,qY 



After calculation we obtain 



H = p 2 {I p + J p ) + q 2 (I q + J q ) + (I p + Iq){J p + J q ) + (a p a q /3 p f3q + a p a q (3 p l3 q ) 

= P 2 (Ip + J p ) + q 2 {I q + Jq) + {Ip + I ? )(Jp + J q ) + 2(lpIqJ p Jq) 1/2 COS^o), 
With 1pQ = 9q-6 p + Lp p - Ifg. 

The Hamiltonian system associated to if is defined on the phase space T 4 xM 4 3 (6 P , 9 g , ip p , <p q \ I p , I q , J p , J q 
by 



(3.1) 



dH 

w 

dH 

&7? 



J, 



dH 

"W 

dH 

dtpj' 



3 =P, 



J = Pi 



Since the Hamiltonian H only depends on one angle (0o), the system (3.1) is completely integrable 
(this is also a consequence of the invariance properties recalled in (1.3)-(1.5)). 

Lemma 3.1. — The system (3.1) is completely integrable. Moreover, the change of variables 



(3.2) 



Iq + Ip, K 2 



Jq + J„, K 3 = I q + J q , K = I q 



4>l = p , V>2 = Pp, ?p3 = (fq- <fp, V'o 

is symplectic: dl A d6 + dJ A dip = dK A dip. 



+ ipp-tp q 



BEATING EFFECTS IN CUBIC SCHRODINGER SYSTEMS 



9 



Proof. — It is straightforward to check that 

K\ = I q + I p , K 2 = Jq + J p and K 3 = I q + J q , 
are constants of motion. Furthermore we verify 

{K 1 ,H} = {K 2 ,H} = {K 3 ,H}=0, 

as well as 

{K 1 ,K 2 } = {K 2 ,K 3 } = {K 3 ,K 1 }=0. 

Moreover the previous quantities are independent. So H admits four integrals of motions that are 
independent and in involution and thus H is completely integrable. □ 

In the new coordinates, the Hamiltonian H reads 

(3.3) H = H(ip ,K ,K 1 ,K 2 ,K 3 ) 

= p 2 (K x + K 2 - K 3 ) + q 2 K 3 + K X K 2 + 2[K{K 3 - K){K 2 -K 3 + K){K X - K)] 1/2 cos Vo- 
We set K\ = K 2 = K 3 = e 2 , and we denote by 

iWo, K ) := H(<po, K , e 2 ,e 2 , e 2 ) = e 2 {p 2 + q 2 ) + e 4 + 2K (e 2 - K ) cos Vo- 
The evolution of (ipo,Ko) is given by 



dH 



dK dip 
Then, we make the change of unknown 

(3.4) Vo(i) = ^(e 2 t) and K (t) = e 2 K(e 2 t). 

An elementary computation shows that the evolution of (tp, K) is given by 



(3.5) 



where 



^ = -2(1- 2K) cosip -- 
K = - 2K(1 - K)sinip 



dK 

dip ' 



H+ = H*(ip,K) = 2K(1 — K) cosip. 



The dynamical system (3.5) is a pendulum whose phase portrait is drawn in Figure 1 and we easily 
deduce that 



Lemma 3.2. - Let 7 > arbitrary small, then the dynamical system ( |3.5| ) admits a periodic orbit 
r 7 := {(^ 7 (i),if 7 (t)) I t e M} of period 2T 7 satisfying (V» 7 (0),if 7 (0)) = (0,7) and (^ 7 (T 7 ), # 7 (T 7 )) = 
(0,1-7)- 
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1 




-tt/2 







tt/2 



Figure 1. Level sets of the Hamiltonian H+. 



Remark 3.3. — For the system (1.7), we can perform a similar analysis, and we get the reduced 
Hamiltonian 

H* = H+W, K) = \{q 2 - P 2 )K + 2K(l - K) cos j> . 
This Hamiltonian is of pendulum type for e arbitrary small iff q = ±p, and with this choice only we 



are able to prove the same result as for (1.1) 



4. Proof of Theorem II. II 



Consider the Hamiltonian H given by (2.4), which is a function of (atj, ctj, (3j, We want to 

prove that, for a good choice of initial datum, the solution of the Hamiltonian system governed by H 
remains close to the solution of the reduced system gove rned by H (or H+ 
We make the linear change of variables given by Lemma 



3.1 



Then H induces the system 



(4.1) 



dH 




f dH 


dH 


9Kj 




ioik - a - , 
da k 


da k 


dH 


< 


■R 911 








I df3 k 





H + Rl + Z42 + Z44 + i?g, 



Denote A = {p, q} and £ = Z\ {p, q}. By Propositions 2.1 and 2.2 we have 

H 

where 
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the polynomial Z\^ contains all fourth order monomials with 2 indices outside A, and Z44 contains 
all fourth order monomials with 4 indices outside A. More precisely 

Z 4 ,2 = {Ip + I q )^Jj + {Jp + J q )^I j + 

+{a p P p + a q f3 q ) ^2 + (®pPp + «?Az) E 
Notice that Ri vanishes when Ik = for all k ^ A. 



Observe that the Kfs aren't constants of motion of (4.1). However, they are almost preserved, and 
this is the result of the next lemma. 

Lemma 1^.1. — Assume that 

(4.2) aj(0),/3j(0) = 0(e), Vj G .4 and aj (0),^{0) = 0(e 2 ), Vj G C. 
Then for all < t < Ce~ 3 , 

(4.3) Ij(t), Jj{t) = 0(e 4 ) w/ien je£, 
and 

(4.4) ^(t) = Ki(0) + O(e 6 )t 

(4.5) # 2 (*) = K 2 (O) + 0(e 6 )t 

(4.6) K 3 (t) = K 3 (0) + 0(e 6 )t, 

where the Kj 7 s are defined by ( |3.2[ ) . 

Proof. — We first remark that by the preservation of the I? norm in each NLS equation, we have 
Ij (t) = 1 i (0) and Y J J (*) = zZ 3 i (°) for a11 * G M ' 



and therefore by using (4.2) 

I n (t) = 0(e 2 ), J n (t) = 0(e 2 ) for all n G Z and for all f G M. 



On the other hand by Propositions |2.1| and |2.2[ we have for n£Z 

(4.7) i n = {l n ,H} = {l n ,Z 4 } + {l n ,R 6 }, 

and the same for J n . 

• To prove (4.3), we compute 

(4.8) { I n , Z 4 } = -i{a n (3 n S - a n (3 n S), 

and { J n , Z4 } = i(a n (3 n S — a n (3 n S). Then by ( |4.7[ ), if we denote by L n = I n + J n we get 

L n = { L n , Rq }. 
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Furthermore, for n ^ p,q all the monomials appearing in { L n ,R,Q } are of order 6 and contains at 
least one mode in C Therefore as soon as (4.3) remains valid, we have L n (t) = C(e 2+5 ) and thus 
|L n (t)| = 0(e A ) + t 0(e ). We then conclude by a classical bootstrap argument that (4.3) holds true 
for t < Ce- 3 . 

• It remains to prove (|4.4l)-(|4.6|). We denote by 



Z% = {I p + I q )(J P + J q ) + (a p a q (3 p (3 q + a p a q f3 p /3, 



qj i 



the fourth order part of the model Hamiltonian. From (4.8), we get 

{ I p , Z%\ = -i(a p f3 p aq(3 q - a p f3 p a q (3 q ) 

= — { Iq: -^4 } = — { Jp, Z%} = { J q , Z%} , 

thus {Kx,Zl} = {lp + I q ,Zl} =0. Similarly, { K 2 , Z% } = and { K 3 , Z% } = 0. Therefore, by 
using fl47l ) we deduce that for all j E {l, 2, 3} 

(4.9) K j = {K h Z A>2 } + {K j ,R & }. 

Then we use that each monomial of Z^^ contains at least two terms with indices j G C Therefore, as 
soon as ((42]) holds, \{Kj, Z 4 _ 2 }| < Ce 6 . Furthermore \ { Kj, R 6 }\ < Ce 6 . Therefore, by ( |49| ), 



K ] {t) = K0) + tO{e & ). 

From now, we fix the initial conditions 

K 1 (0) = e 2 , K 2 (0) = e 2 , K 3 (0) = e 2 , 
and |oy(0)|, 1^(0)1, |^(0)|, |^(0)| < Ce 2 for j^p,g. 



□ 



(4.10) 



Let H be given by (2.4). Then according to the result of Lemma 4.1 which says that for a suitable 



long time we remain close to the regime of Section |3j we hope that we can write H = Hq + R, where 
R is an error term which remains small for times < t < e -3 . 

We focus on the motion of (ipo, Kq) and, as in the previous section, we make the change of unknown 

(4.11) Tp (t) = i)(e 2 t) and K (t) = e 2 K(e 2 t), 

and we work with the scaled time variable r = e 2 t. Then we can state 



Proposition — Consider the solution of (4.1) with the initial conditions (4.10). Then (i{),K) 
defined by (4.11) satisfies for < r < e^ 1 



(4.12) 

where H+ is the Hamiltonian 



- - 9 lt + 0(e 2 ) 
2K(1 -K) cosip. 
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Proof. — First recall that H = H{ipQ, Kq, K\, K2, K3) is the reduced Hamiltonian given by (3.3). By 
Propositions 2.1 and 2.2 we have 

(4.13) H = H + R I + Z 4:2 + Z 4A + R 6 . 
Thanks to the Taylor formula there is Q so that 

H(^ ,K ,K 1 ,K 2 ,K 3 ) = H(^K ,e 2 ,e 2 ,e 2 ) + Q 

(4.14) = H + Q. 
Thus, by $~Tty and (J4T4]) we have H = H + R with 

R = Q + Ri + Z42 + Z4A + i?6- 



By (4.1), (ipo,Ko) satisfies the system 

Mt)= -§f (Mt),K (t),...) 

Ko(t)= §^ (Mt),K (t),...), 

where the dots stand for the dependance of the Hamiltonian on the other coordinates. Then, after 
the change of variables (4.11) we obtain 

^(r) = 



-iff(^(r),^(r),...) 



Now write H = Hq + R and observe that Ho(ip,e 2 K) = C e + e^H^tp, K). As a consequence, (ip,K) 
satisfies 

dH* 1 dR(iP,e 2 K,...) 



8K 
dHi, 



dK 



1 dR(-4>,e 2 K,...) 
dtp + dtp 

Thus it remains to estimate d, l i ) R(ip,e 2 K, . . .) and djcR^ip, e 2 K, . . . ). Remark that ^ and -ff are 
dimensionless variables. Thus, if P is a polynomial involving k internal modes, (ay, ay, /3j, (3j)j^, and 
^ external modes, (ay, ay, f3j)j e c, we have by using Lemma 4.1 

d 1 pP(i;,e 2 K, ...) = 0(e k+2t ), d K P^,e 2 K, ...) = 0(e k+2t ). 

As Ri contains only monomials involving at least one external action (Ik)k<£A we S e ^ 

d^R T (<iP,e 2 K,...) =0(e 6 ), d K Ri{$, e 2 K, . . . ) = 0(e 6 ), 

d<pZ 4i2 (^,s 2 K,...) =0(e 6 ), d K Z A ^,e 2 K,...) = 0(e 6 ), 

a^z^^K,...) =o(e 6 ), a K z^,e 2 K,...) =o(e 6 ), 

c^ 6 (V,£ 2 K,...) =0(e & ), d K R 6 (^,e 2 K,...) =0(e 6 ). 

On the other hand, by construction Q reads P\/S.K\ + P2AK2 + P3AK3 where Pi, P2 and P3 are 
polynomials of order 1 in Kq, K\, K~2, K3 and e 2 while AKj denotes the variation of Kf AKj = 
Kj — Kj(0). Using again Lemma 4.1 we check that for < r < e _1 

fyQ = 0(e 2 ), d K Q = 0(e 2 ), 

hence the result. □ 
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We now consider the solution (tpy,Ky) of (3.5), described in Lemma 3.2, which is issued from the 
initial condition {ip y, i£y )(0) = (0,7) for some 7 such that e 1 / 2 < 7 <C 1 and we compare it with the 
solution (ip,K) of (4.12) issued from the same initial datum: 



For all < r < e 1 we have 

0;,K)(T) = (^,K J )(T) + O(e 2 )T. 



Lemma 4-3. 

(4.15) 

Proof. — Consider the system ( |3.5| ) and the open domain U = {—it, n) X (0, 1). By the Arnold Theorem 
(cf. [1, p. 113], see also [HJ Lemma 4.3]), this Hamiltonian system admits action-angle coordinates 
(L, a) = Q(ip,K) defined onW\ {(0,0)} by a C 1 symplectic map <J? satisfying that uniformly on any 
compact U CU\ {(0,0)}: 



|d$|| < C, 



|d$ 



< C . 



Then we obtain that for < r < e 1 

d 



dr 



(L,a) 



dr 



${i>,K) = 6$($,K).{j),k) 



dip ' dK 



+ 0(e 2 



Therefore there exists L* € M so that L{t) = L± + C(e 2 )r and if we define uj* = — ^^(L*), we obtain 
a(r) = + 0(s 2 )t. Notice that by construction $(^ 7 (r), if 7 (r)) = (L*,ovr) for all r£i Next, 
as d<3> _1 is bounded, we get 

^,K)(t) = $- 1 (L(r),a(r)) = S" 1 (L*, ^r) + 0(e 2 )r 

= (^ 7 ,^ 7 )(r) + 0(e 2 )r. 



With the choice e 1 / 2 < 7, the remainder term in (4.15) is so that |0(e 2 )r| <C 7 < K-y for t < e 1 . 



□ 



Proof of Theorem \l.l\ — As a consequence of Lemma 4.3, the solution of (4.1), with initial datum 
Q4.10P and (ip , K )(0) = (0,e 2 j), satisfies for < t < 

K (t) = e 2 K 1 (e 2 t) + 0{eH) 

Mt) = *Py(e 2 t) + 0(eH), 

and with the condition e 1 / 2 < 7 we obtain (1.2). 

We now compute the period 2T 7 . From the expression of the Hamiltonian H+, we infer 



K = 2K(1 - K), 



h? 



[2K{1-K)) 



2 > 



where h = 27(1 — 7) > 0. Thanks to the symmetries of H+, T 7 /2 is the travel time for the solution 
between (0,7) and (— cos _1 (2/i), 1/2). In the interval [0,T 7 /2], the function K is strictly increasing, 
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hence invertible, so we can write the time t as a function of K, and this implies that 

(4.16) T 1 = 2f l — ^ K 

7 7 7 y/(2K(l - K)) 2 - h 2 

Next, we estimate T~. It is easy to check that there exists C > such that for all 7 < K < i, we have 



K 2 (l 



^)2_ 7 2 (1 _ 7) 2> C( ^2_ 72)> 

Hence by integration in (4.16) we deduce that there exists C > so that 
(4.17) T 7 <-C In 7. 



□ 



5. Proof of Theorem 11.31 



Fix a > 1 and consider the system (1.1). We fix the set of internal modes: p = 0, q 3> 1 and 



e := e 2* 



Then we consider the first equation in (1.1) as a linear time-dependent Schrodinger equation 

id t u + dlu + V q {t,x)u = 0, (t,x)£RxE 1 , 
with potential V q = — e~ ql/a \v\ 2 . 

The regularity properties of the solutions given in Theorem [lT] imply that, on [0, S 1 the function 

V q is smooth in time and real analytic in space. In order to construct the sequence of initial conditions 
announced in Theore m |1.3[ we have to ensure uniform bounds w.r.t. the integer q E N in the Gevrey 
class ^(S 1 ), given in ( |1.9[ ). We have 

V g (t,x) = e-« lA > (i) + v q (t)e iqx + e-^ 1/a r v (t,x)\ 2 , 

where r v (t, .) is an analytic function, whose norm is uniformly bounded with respect to \t\ < and 
q G N. We then compute the Fourier coefficients Vj of V g (t, .). 

The dominant coefficients are labelled by the indices 0, q and — q: for them, we have 

\Vo\<e-« 1/a (\v \ 2 + \v q \ 2 + e-y i/a R), 

\V q \<e-* 1/a (\v \\v q \+e-y i/a R'), 

\V- q \<e-« 1/a (\vo\\v q \ + e-y i/a R"), 

where R,R',R" are uniformly bounded w.r.t. q G N and < t < e~ 3 . Since vq and v q stay (in 
modulus) between and 1, estimate (1.9) is obtained for these indices. 

The coefficients Vj for j ^ —q,0,q decay much faster in general: using the analyticity of r v (t, .), 
and the fact that for every q G N and t G Z 

\\ l/a + p\t-q\ >c\£\ l/a , 



4 

we obtain 

|VSI < Ce-^'l 1/a . 
Once again, this estimate is uniform w.r.t. q G N and < t < e -3 . 
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Choose initial c ond itions so that (tp(0), K(0)) = (0,7) for some < 7 <C 1. In order to apply the 
result of Theorem 1.1 we must have T 7 < e . Therefore we impose 



T — — I < 

q -~ e 2 e 3 ' 

which leads to e~z q 1 < smce ^7 — C| m 7l- We fix the IP-norm of the initial condition with 

the choice 7 = q~ 2s (observe that for q 3> 1, e < 7 2 , so that we are in the conditions of application 

l„l/a 



1.1 



then the previous constraint becomes e 2' < 2s\nq ^ bence is satisfied for q large 



of Theorem 

enough. From (4.17), we get 

T q < C|ln 7 |e- Z = 2sCe r '"lng. 
Now, the growth rate of ||u||fp between t = and t = T q is bounded from below by C'7~ 1//2 , where 



C is independent of s. So we have, by (1.6), 



\H0)\\h° - (l + 2a 8 ) aaK q 



as 



Note that the constant C ajS := n+l^pi S oes to as s or a goes to infinity. □ 

Remark 5.1. — If we choose a different e, as for example e = exp (-(lng) 1+K ), with k > 0, wc have 
that for all s > 

VgeN, Vt€[0,T,], ||V^(t, < C s ,a , 

and we obtain the growth 

\\u( T q)\\H* 
\\u{0)\\ H s 

that is, a sub-polynomial growth. 



> Cexp (aClnT,) 1 /^*)) , 
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